Quantum optical weak measurements can visualize photon dynamics in real time 
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An experiment is proposed to visualize stroboscopically in real time the dynamics of a photon 
oscillating between two cavities. The visualization is implemented by a sequence of weak measure- 
ments (POVM), which are carried out by probing one of the cavities with a Rydberg atom and 
detecting a resulting phase shift by Ramsey interferometry. This way to measure the number of 
photons in a cavity was experimentally realized by Brune et al. . We suggest a feedback mechanism 
which minimizes the disturbance due to the measurement and enables a detection of the original 
evolution of the radiation field. 

PACS numbers: 42.50.-p, 03.65.Ta, 32.80.-t, 03.67.-a 



There is much experimental progress in trapping single 
atoms, ions and photons. It has recently been reported 
that it is possible to trap individual atoms with a sin- 
gle photon in a cavity and to reconstruct the trajectory 
of the atom Q. We want to contribute to this rapidly 
growing field of quantum visualization by proposing an 
experimental realization in a very "clean" setting. In 
a preceeding paper Q we have shown theoretically and 
numerically that it is feasible to monitor in real time a 
dynamical process occurring in an individual two-level 
system with state vector 



\i>{t)) =Ci(t)|¥>i) +c 2 (t)\tp 2 ) 



(1) 



Our aim in the following is to describe an experimen- 
tal set up for quantum visualization which serves to reg- 
ister the time behavior of |c2(i)| 2 approximately, while 
only weakly influencing the original dynamics of \ip(t)). 
Here we assume the original dynamics to be known be- 
fore the measurements. Why should we want to measure 
a theoretically known dynamics? Because we want to 
verify that our measurement procedure works. It is a 
non-trivial task to measure | £2 (i) | 2 _ a quantity which 
refers to an ensemble - by a one run measurement on a 
single system in real-time. In order to achieve this goal 
we need only little a priori information about the original 
dynamics. In the concrete example of Rabi oscillations 
presented in this paper we succeed with the proposed 
measurement scheme knowing only the order of magni- 
tude of the oscillation period. Our long-term objective is 
to devise a measurement scheme for the real time moni- 
toring of partly unknown dynamics. 

In classical physics it is possible to track the evolu- 
tion of an individual system without disturbing it sub- 
stantially. How can this aim at least approximately be 
achieved for quantum systems? Since projection mea- 
surements, which are also called sharp measurements, 
severely alter the original motion of \ip(t)), they are not 
suitable in an one shot situation where only a single real- 
ization of this motion is available. An exception are QND 
schemes. They have the disadvantage to require for non- 
trivial dynamics an observable with a continuous spec- 
trum H, which does not exist for two-level systems. In 
case of a two-level system one needs instead weak (or un- 



sharp) measurements (POVM measurements) by which 
the state of the system is less disturbed but nevertheless 
some information about the state is provided. A single 
weak measurement can be realized by suitably entangling 
the two-level system with a quantum meter via a unitary 
transformation (premeasurement) followed by a projec- 
tion measurement on the meter. The latter supplies a 
measurement result, which is read off. 

To track the development of | £2 (*) | 2 in time, a sequence 
of weak measurements is necessary. The corresponding 
series of measurement results can then be appropriately 
processed in real time to give the final measurement read- 
out. Two conditions may thus be fulfilled simultaneously: 
i.) The back action of the measurements does only mod- 
erately disturb the original dynamics of the system given 
by the evolution of \ip(t)) and ii.) the variance of the mea- 
surement results is small enough to enable a reliable esti- 
mate of the original time behavior of | £2 ) | 2 - Of course 
there is no information gain without disturbance of the 
state. In fact the greater the information gain the greater 
is the change of the state due to the measurement. 

In this paper we sketch an experiment to visualize Rabi 
oscillations with frequency When fulfilling the con- 
ditions i.) and ii.) we have some freedom in adjusting 
the measuring apparatus. For example Sli? does not have 
to be known exactly beforehand. An apparatus tuned to 
any frequency out of the interval [0.75 flfj , 1.5 Qr] would 
reveal SIr as the actual frequency (see below). 

In cases in which a sequential measurement may be ap- 
proximately treated as a continuous measurement it can 
in the selective regime be described by a stochastic mas- 
ter equation . This powerful calculational tool enables 
an optimal exploitation of the information contained in 
the readout of the continuous measurement || . 

We are considering instead a series of well separated 
weak measurements. Such a series was employed to carry 
out a QND measurement of small photon numbers in an 
experiment of Brune, Haroche et al. ||, which was the- 
oretically analyzed in Q and experimentally realized in 
||. While the weakness of the measurements has been 
considered as an obstacle there, it turns out to be an ad- 
vantage when it comes to the detection of dynamics. The 
experimental setup we sketch in the following is based on 
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the Brune-Haroche experiment. We have added a second 
cavity and a feedback mechanism. The system consists of 
one photon with frequency u> shared by two equally con- 
structed, coupled cavities C\ and C 2 . One could think of 
two identical cavities connected by a waveg uide (cf. (|) 
or a transmissive mirror. For calculations the cavities are 
assumed to have infinite damping time. A justification 
will be given bellow. Their coupling is modeled by the 
interaction Hamiltonian of the Jaynes-Cummings type: 



H = %g{aia\ + a\a 2 ) 



(2) 



with coupling constant g. In the interaction picture, 
which we are going to use, (§) is the full Hamiltonian. 
The indices refer to the cavity numbers. Such a coupling 
of two cavities has also been considered in 0. The 
photon which is delocalized over the two cavities can be 
described as a superposition of two states: 



|^(i)} = ci(i)|l,0)+c 2 (t)|0,l). 



(3) 



The first and the second slot in the ket represent the 
number of photons in cavity C\ and cavity C 2 respec- 
tively. For the initial state \i/j(t = 0)) = |1,0) we find 
Rabi-oscillations with the Rabi-frequency £Ir := 2g 



\c 2 (t)\ 2 = sm 2 {gt) . 



(4) 
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FIG. 1. Relevant levels of Rydberg atom 

It is our goal to measure this original evolution of 
\c 2 (t)\ 2 in real time by probing the coupled cavities with 
atoms. We sent a Rydbergatom with three effective en- 
ergy levels g, e, i (see Fig.[jj) and velocity v through the 
first cavity C\ (cp. |j]|7)). The passage time L c /v ( L c is 
the cavity length) is assumed to be much shorter than the 
period Tr := 2tt/Q,r = n '/ g of the oscillations of | £2 (i) | 2 ■ 
Then the coupling of the two cavities is negligible during 
the time the atom spends in the cavity. The detuning of 
the atomic transitions with respect to the frequency of 
the cavity mode to is such that the interaction between 
the atom and C\ is dispersive and only the energy levels 
e and i suffer an appreciable dynamical stark shift. Pro- 
vided the atom enters the cavity in in a superposition of 



states \g) and |e), the effective Hamiltonian reads (cp. 
eqn. (16) in [§): 



nn 2 . . . . t 
—jr- |e)(e| ® a[ai , 



(5) 



where S := uj — uOi e and f2 = f2(r) is the Rabi frequency 
averaged over the path of the atom through the cavity. 
With (||) the state of the enlarged system composed of 
the atom and the photon field changes according to 



(c e |e) + c g \g)) ® 



\e)®U Cl \il>) 
c g \g) ® 



(6) 



with C/ci being diagonal in the basis |1, 0) and 0, 1): 

U Cl :=e- l£l |l,0)(l,0| + |0,l)(0,l| , (7) 

and £1 :— t- — ■ \ip) represents the state of the 1-photon- 
field probed by atoms. The net effect of the atom-field 
coupling described by the interaction Hamiltonian (j|) is 
that only the amplitude of the alternative |e) <£> |1, 0) suf- 
fers a phase shift e~ l£1 while the amplitudes of the other 
quantum alternatives remain unchanged. 
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FIG. 2. Experimental setup 

Phase shifts between several quantum alternatives may 
be measured by interferometry. As proposed in || it is 
convenient to use the Ramsey method of separated os- 
cillatory fields (see Fig. ||). To this end a Rydberg atom 
is initially prepared in state \g). Before entering cavity 
C\ the state of the atom is transformed into a superpo- 
sition of states |e) and \g) by entering a cavity which 
contains a first classical oscillatory microwave field Ri 
with frequency u> r . In the cavity C\ the atomic state 
becomes entangled with the state of the cavities as dis- 
cussed above. After leaving C\ the atom crosses a second 
cavity with a classical microwave field R 2 which is in 
phase with R\ and positioned at the distance L from it. 
The total state change of duration St amounts to (cf . ) 
l*(*o)) - * |^(*o + St)), where the product state before 
the atom enters C\ is given by 

l*(*o)> = \g) ® |^(*o)> = \g) ® (ci(t )|i,0) + ca(t o )|0, l)) 

(8) 

and the final entangled state equals 
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\9(t + Sr)) = |e>® (ttfd(io)|l > 0) + «5c 2 (t )|0,l» 
+ b> «> (ufci^o)!!, 0) + ufc 2 (to)|0, 1» 



The coefficients in (|^) are given by 
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with e = characterizes the Ramsey fields and 

depends on the length l r of each of the corresponding- 
cavities and the effective Rabi-frequency Q r inside these 
cavities: vq :— 2l r Q r /it. <po := (uj r —uj eg )^ is the phase 
shift which is induced by the Ramsey cavities in the case 
v = Vq. An analogous result was obtained in eqn. (A7) of 
JjJ for the initial atomic state being \e). Eqn. (||) shows 
that the meter states |e) and \g) couple in general to both 
cavity states |1, 0) and |0, 1). This is a characteristic trait 
of a weak measurement. 

After the atom has left the second Ramsey field R 2 its 
energy is finally detected in a projection measurement 
by field ionization counters D e and D g . The state of 
the composite system after a measurement with result 
I e {e,g} reads |*j(* + St)) = \l) <g) \ipi(t + St)} with 
photon state 



Mto + ST)) = |ui|ci(t )|l,0) 



(11) 
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and Uj — \u l j\e lx i for j G {1,2}. Here a global phase 
factor has been omitted. The probability to obtain the 
related measurement result I is given by the expecta- 
tion value of the corresponding projector: prob(7) = 
(( \l)(l\ ® 1) )*( to +i T ). Eqn. ( [il| ) shows that after the 
measurement the photon is in general not localized in one 
of the cavities. The disturbance of the photon state due 
to the measurement may be small. Because of the Rabi- 
cvolution between the measurements this set up repre- 
sents no QND measurement of the photon number as it 
has been in the Brune-Haroche experiment. 

Referring to the photon field only, the change of its 
state due to a single measurement with result I can be 
expressed by an operation Mf. | "0 (*o ) ) —* \tf>l{to + St)) = 
Mi\tp(to)). Like all bounded operators, Mi can be written 
as "phase" times "modulus" (polar decomposition) 

Mi = Ui\Mi\ (12) 

with unitary transformation 

Ui := |l,0}(l,0|+e l ^-xi)| 05 i)( 0jl | (13) 



and positive operator 

\M t \ := lull 11,0X1,01 + 14110,1X0,11 . 
The probability to obtain the outcome I is then: 



(14) 



prob(Z) - (MjM l ) v , (to) = < \M,\ 2 )^ {to) . (15) 

The Ei := \Mi\ 2 is also called effect. In this way we 
obtain e.g. for the probability to measure the energy e: 



Pe 



Pi|ci| 2 +p 2 \c 2 \ , where p 3 



,e|2 



is fixed by (10). 



The effects have the property E e + E g — 1 and gener- 
ate a positive operator valued measure (POVM). In the 
special case where u\ — u 9 2 = 1 and 



g 



0, the 

operation Mi = Ei is a projector. If on the other hand 
Ei = 1, no measurement at all has taken place but only 
an unitary development (Mi — Ui). These two cases are 
the extremes of a sharp and a totally unsharp measure- 
ment. By varying the parameters v , vq , tpo and e of the 
setup all degrees of "weakness" between these two ex- 
treme cases as well as the extremes themselves can be 
reached. 



Eqn. (|15[) shows that the information obtained by the 
generalized measurement is solely contained in |M/|. This 
part of the operation Mi in (12) represents at the same 
time the unavoidable minimal disturbance of the system 
by the measurement. But our set up causes in addition 
by means of Ui a purely unitary or Hamiltonian evolution 
of the state, which modifies the photon state without be- 
ing necessary for the extraction of information. Since we 
want to disturb the original state motion as little as possi- 
ble, we have to install a Hamiltonian feedback mechanism 
which compensates Ui given by (|l3|). Such a procedure 
has already been proposed by Wiseman jll| . In our case 
feedback can be implemented by modifying the set up as 
follows: 

After a measurement beginning at an arbitrary time 
t = to with outcome I an atom prepared in state |e) is 
sent through the second cavity C 2 - As in the case where 
an atom crosses cavity C\ the unitary evolution is again 
governed by the dynamical Stark effect, with the only dif- 
ference that now the energy shift depends on the number 
of photons in C 2 instead of C\ : 



\e) ® \ipi(t + St)) 



®Uc 2 \Mto+ST)) 



with 



U C2 := |l,0)(l,0| + e- l£2 |0,l)(0,l| 



(16) 



(17) 



The combined influence of the measurement and feedback 
leads to 



Uc 2 Mto + ST)) = UcMMi\mto)) 



(18) 



The condition for compensation of Ui in di~3| ) is therefore 
Uc 2 Ui = 1 <^> £2 = X l 2 ~ Xi: where Xj ma y be obtained 
from (|l^). This condition demands that the compensat- 
ing phase e 2 — jj^j (f denotes the feedback) has to be 
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chosen depending on the measurement outcome I. We 
see two ways to vary 62- One is to select an appropriate 
velocity Vf of the feedback atom sent through the upper 
cavity Ci. The other possibility consists in setting up a 
suitable detuning 5f. This can be done by shifting the 
atomic energies by means of an static electric held in the 
cavity C2 cp. pig] . Please note that it makes no difference 
whether the atom sent through cavity C2 is thereafter 
measured or not because the composite system after the 
interaction is in a product state. 

In order to reach our final aim of monitoring the orig- 
inal Rabi-oscillations of | £2 (t) | 2 , a sequence of measure- 
ments at times t n = nr with n = 1,2,3... has to be 
carried out. Between two consecutive measurements the 
system evolves undisturbed according to the Hamiltonian 
(||) . The resulting total evolution of the system is given 
by C2(t) instead of C2(t). To process the data obtained 
in the single measurements we first of all divide the se- 
quence of results with values e and g into groups of N. 
From each so called "AT-series" we extract the relative 
frequency r :— N e /N of the number N e of e-results. It 
turns out that its expectation value £(r) is related 
to the value which | C2 1 2 assumes immediately before the 
start of the iV-series by 



C2 



£(r) - Pi 
Ap 



(19) 

of @. In a se- 



with Ap := P2 — P\ = 1 2 — I 
quence of measurements on a single radiation field we 
do not have access to the expectation value £(r), which 
refers to an ensemble. Instead we insert r(to) into the 
right hand side of equation (|l9|) and obtain thus a "best 
guess" of |c2 (^o) | 2 at time to when the first measurement 
of the respective N-series began: 



G a (fo) 



r(tp) - pi 
Ap 



(20) 



The possible values of G2(fo) are distributed around 
C2(^o)| 2 aim mav be negative. This estimation of |c2(io)! 2 
can be good only if the duration of the N-series Nt is 
much smaller than the period Tr of the oscillations of 
the system. The sequence of G2 at various times serves 
as the final readout of the sequential measurement. 

We have two competing influences on the system: The 
strength of the original dynamics is proportional to g or 
fin = 2tt/Tr . The measurements on the other hand 
hinder this dynamics the more the stronger they are and 
the quicker they are repeated with frequency 1/r. The 
Zcno effect demonstrates this! A measure for the distur- 
bance due to the sequence of measurements is the deco- 
herence time To, the time after which the off-diagonal 
elements of the density matrix (in the \tfi), \f2) basis) 
have decayed to 1/e of their original value. In a sub- 
sequent paper we will show that To = 8 P °^~^°V with 

po := (pi+p2)/2. In order to have a high resolution in the 
sequential measurement it would be desirable to have a 



small decoherence time (To is proportional to the time it 
takes to distinguish between the states |y>i) and \ip2)i C P- 
0] ) . On the other hand the influence of the measurement 
should not dominate the evolution. We therefore require 
the decoherence time to be maximally as great as the 
Rabi time. A numerical analysis showed that for our 
purpose a favorable balance of information gain and dis 
turbance is obtained if the so called fuzziness / 
is adjusted to be close to one: / 
to choose / in the interval 0.75 < / < 1.5. The exper- 
imental parameters e, ipo, vq, v and r have to be fixed 
correspondingly. 



2Tr 

1. In fact it suffices 
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FIG. 3. Under a sequence of appropriate weak measure- 
ments the measurement readout G2 (grey curve) is corre- 
lated with the state evolution | C2 ) | 2 (black curve). This 
becomes evident after noise reduction (dashed curve) of 
the readout G2, which was carried out taking into ac- 
count approximately 12 Rabi-cycles. Parameter values: 
^ e [1.25 - 0.125,1.25 + 0.125], e = 0.068 • w, (p = tt, 
r = 0.002T/J and A^ = 25, which lead to an average fuzzi- 
ness of / = 0.98. 

We have simulated numerically all the processes de- 
scribed above including the feedback. In a realistic ex- 
periment the velocity v of the probing atoms and the 
feedback atoms will vary from one single measurement 
to the other. We have accordingly tolerated the veloci- 
ties v to fluctuate uniformly by ±10% about the desired 
mean value. The resulting dynamics of the state under 
the influence of stroboscopically applied weak measure- 
ments is given by the |c2(i)| 2 -curve (black) in Fig. ||. The 
measurement readout which is defined as best guess G2 (t) 
of I C2 (i) 1 2 (grey curve) has been further processed to the 
noise reduced G2-curve (dashed). The noise reduction 
procedure consists essentially in a time-averaging of the 
readout. Details are described in Appendix C of 0. 

We find a high correlation including the phase between 
the noise reduced G2-curve and the |c2| 2 -curve. The ac- 
tual evolution of the state is therefore well monitored in 
time. The |c2| 2 -curve reflects the fact that the original 
Rabi-oscillations have been disturbed by the measure- 
ment, though they are only slightly modified. 

Fig. ^ shows the power spectrum of the |c2| 2 -curve 
(black) and the measurement readout G2 (grey). Both 
curves are peaked at the Rabi frequency £Ir. 
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FIG. 4. Power spectra of |c 2 (t)| 2 (black) and G 2 (t) (grey). 

All parameters used in Fig. || and || are realistic (cf. 
0) apart from the number of Rabi cycles and the effi- 
ciencies of the detector and the feedback mechanism. In 
what follows, we are going to address these remaining 
problems. 

The number of Rabi cycles which can be monitored 
depends crucially on the mean lifetime tcav of the cou- 
pled cavities which has so far been assumed to be infi- 
nite. At the required frequency of 50 GHz a icav = 0.1 
s seems to be feasible (cf. 0). The time of flight of 
an atom with velocity v — 1000 m/s through cavity C\ 
amounts to Lc/v = lOfis. We therefore assume that a 
time t — 100/is between two consecutive measurements 
suffices to send a feedback atom through cavity Ci be- 
tween the measurements. Then 1000 measurements can 
be made within the mean lifetime icav of the cavities. 
For Tr/t = 500 as assumed in Fig. |], this corresponds 
to monitoring the first two Rabi cycles. As can be seen in 
Fig. |^ the readout of the measurement G2 (grey curve) os- 
cillates about the evolution of the component | C2 1 2 (black 
curve), so that also for two Rabi cycles the evolution of 
I C2 1 2 can be approximately recovered by averaging G2 
over an appropriate timescale. 

A more serious problem arises when some of the feed- 
back atoms are not sent through cavity C2 and thus the 
unitary part of the back-action Ui is not always compen- 
sated. In our simulations we obtained still good results 
with 4 percent of the feedback atoms missing while for 
more than 7 percent the evolution of | C2 1 2 is significantly 
disturbed. Since the preparation of the Rydberg atoms 
in the Haroche experiment is a random process, feedback 
atoms might be missing. This problem could be solved by 
sending a high number of atoms at once through cavity 
C2, each with large detuning, i.e., causing a small phase 
shift of the radiation field. The mean number n of atoms 
sent in one burst should lead to a total phase shift, which 
compensates Ui. Since the standard deviation divided by 
the mean number An/fi can be made arbitrarily small by 
increasing fi, U\ can be compensated very precisely. 

Another crucial point is the detector efficiency. The 
measurement results are not so sensitive to the loss of in- 
formation. It is the disturbance of the Rabi oscillations 
as a consequence of not knowing how to prepare the feed- 
back atoms which mainly hinders the visualization of the 
photon dynamics. The problem can be solved by em- 



ploying a detector with more than 96% efficiency (see 
above). If only detectors with moderate efficiencies are 
available there is another solution. It consists in choos- 
ing the experimental parameters in (^) such that the 
measurement statistics become Poissonian: 



0, u% 



e, u% 



1, and u\ 



(21) 



where e = sui(ttvo/2v) should be small. In this case we 
do not need feedback because Ui = 1 for both I = e 
and I — g . Most measurement results will be "g ", in 
which case the change of the state will be of higher order 
in the small parameter e. The result "e " is likely only 
when I C2 1 2 ~ 1 and then the state change is also very 
small. The disadvantage of the Poissonian method is that 
the quality of the monitoring suffers from the Poissonian 
statistics. The measurement results essentially indicate 
only the parts of the state evolution where | C2 1 2 ~ 1. For 
detector efficiencies higher than 96% the method with 
feedback atoms shows clearly better results. For efficien- 
cies below 96% it is preferable to avoid the necessity of 
feedback atoms. Fig. [| displays a simulation of a se- 
quence of measurements with Poissonian statistics for a 
detector efficiency of 60%. Roughly the same percent- 
age of the maxima of the Rabi oscillations are indicated 
by the readout G2. As in the case with feedback we al- 
lowed the velocities of the Rydberg atoms to fluctuate 
uniformly by ±10% around the mean value v. 




FIG. 5. Simulation of a secfuence of weak measurements 
with Poissonian statistics and a detector efficiency of 60%. 
The back-action of these measurements has no unitary part. 
The dashed and the black curve represent G2 and |c2(t)| 2 re- 
spectively. Parameter values: ~ € [20 — 2, 20 + 2], e — n^, 
0, r = 0.002T/J and TV = 25, which lead to an average 



<P0 = U, T 

fuzziness of / : 



2.04. 



To sum up: Using the feedback mechanism for detector 
efficiencies higher than 96% the original Rabi-oscillations 
are well visualized in phase and frequency. For moder- 
ate detector efficiencies an alternative method leads to a 
moderate visibility of the Rabi-oscillations. 

We wish to thank the referee for his suggestion to use 
groups of atoms for the feedback mechanism in order to 
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improve its accuracy. 
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Konstanz. 
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